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Abstract 

Pseudogap physics in strongly correlated systems is essentially scale dependent. We generalize 
the dynamical-mean field theory (DMFT) by including into the DMFT equations dependence on 
correlation length of pseudogap fluctuations via additional (momentum dependent) self-energy S^. 
This self - energy describes non-local dynamical correlations induced by short-ranged collective 
SDW-like antiferromagnetic spin (or CDW-like charge) fluctuations. At high enough temperatures 
these fluctuations can be viewed as a quenched Gaussian random field with finite correlation length. 
This generalized DMFT+Sk approach is used for the numerical solution of the weakly doped one- 
band Hubbard model with repulsive Coulomb interaction on a square lattice with nearest and next 
nearest neighbour hopping. The effective single impurity problem is solved by numerical renormal- 
ization group (NRG). Both types of strongly correlated metals, namely (i) doped Mott insulator 
and (ii) the case of bandwidth W < U (U — value of local Coulomb interaction) are considered. 
Densities of states, spectral functions and ARPES spectra calculated within DMFT+Ek show a 
pseudogap formation near the Fermi level of the quasiparticle band. We also briefly discuss effects 
of random impurity scattering. Finally we demonstrate the qualitative picture of Fermi surface 
"destruction" due to pseudogap fluctuations and formation of "Fermi arcs" which agrees well with 
ARPES observations. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.27+a, 71.30.+h, 74.72.-h 
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I. INTRODUCTION 



Pseudogap formation in the electronic spectrum of underdoped copper oxides^ is espe- 
cially striking anomaly of the normal state of high temperature superconductors. Despite 
continuing discussions on the nature of the pseudogap, we believe that the preferable "sce- 
nario" for its formation is most likely based on the model of strong scattering of the charge 
carriers by short-ranged antiferromagnetic (AFM, SDW) spin fluctuations^. In momen- 
tum representation this scattering transfers momenta of the order of Q = (-,-) (a — lattice 
constant of two dimensional lattice). This leads to the formation of structures in the one- 
particle spectrum, which are precursors of the changes in the spectra due to long-range 
AFM order (period doubling). As a result we obtain non-Fermi liquid like behavior (dielec- 
trization) of the spectral density in the vicinity of the so called "hot spots" on the Fermi 
surface, appearing at intersections of the Fermi surface with antiferromagnetic Brillouin zone 
boundary 2 . 

Within this spin-fluctuation scenario a simplified model of the pseudogap state was 
studie d 2 ' 4 ^ under the assumption that the scattering by dynamic spin fluctuations can be 
reduced for high enough temperatures to a static Gaussian random field (quenched disor- 
der) of pseudogap fluctuations. These fluctuations are defined by a characteristic scattering 
vector from the vicinity of Q, with a width determined by the inverse correlation length of 
short-range order k = 

Undoped cuprates are antiferromagnetic Mott insulators with U 3> W (U — value of local 
Coulomb interaction, W — bandwidth of non-interacting band), so that correlation effects 
are actually very important. It is thus clear that the electronic properties of underdoped 
(and probably also optimally doped) cuprates are governed by strong electronic correlations 
too, so that these systems are typical strongly correlated metals. Two types of correlated 
metals can be distinguished: (i) the doped Mott insulator and (ii) the bandwidth controlled 
correlated metal W ~ U. 

A state of the art tool to describe such correlated systems is the dynamical mean-field the- 
ory (DMFT) fr7i9i?i 1( ? . The characteristic features of correlated systems within the DMFT are 
the formation of incoherent structures, the so-called Hubbard bands, split by the Coulomb 
interaction U, and a quasiparticle (conduction) band near the Fermi level dynamically gen- 
erated by the local correlation o 6 i 7 i 8 i 9 i 10 . 
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Unfortunately, the DMFT is not useful to the study the "antiferromagnetic" scenario of 
pseudogap formation in strongly correlated metals. This is due to the basic approximation 
of the DMFT, which amounts to the complete neglect of non-local dynamical correlation ef- 
fects. The aim of the present paper is to describe the main results of a semiphenomenological 
approach, formulated by us recently to overcome this difficultyii. 

The paper is organized as follows. In section |H] we present a formulation of the self- 
consistent generalization we call DMFT+Sk which includes short-ranged (non-local) cor- 
relations. Section IIHI describes the construction of the k-dependent self-energy, and some 
computational details are presented in section II V Al Results and a discussion are given in 
the sections HV1 and IV! 



II. INTRODUCING LENGTH SCALE INTO DMFT: DMFT+S k APPROACH 

Basic shortcoming of traditional DMFT approac h^!?!?! 1 ? is the neglect of momentum 
dependence of electron self-energy. This approximation in principle allows for an exact so- 
lution of correlated electron systems (in infinite dimensions) fully preserving the local part 
of the dynamics introduced by electronic correlations. To include non-local effects, while 
remaining within the usual "impurity analogy" of DMFT, we propose the following proce- 
dure. To be definite, let us consider a standard one-band Hubbard model. The extension 
to multi-orbital or multi-band models is straightforward. The major assumption of our 
approach is that the lattice and Matsubara "time" Fourier transform of the single-particle 
Green function can be written as: 

Gk(iuj) = —— — 1 uj = nT(2n+l) } (1) 

iu + n — e[k) — l;{iuj) — ^(iuj) 

where Yi{iu) is the local contribution to self-energy, surviving in the DMFT, while T,k(iuj) 
is some momentum dependent part. We suppose that this last contribution is due to either 
electron interactions with some "additional" collective modes or order parameter fluctua- 
tions, or may be due to similar non-local contributions within the Hubbard model itself. 
To avoid possible confusion we must stress that ^(iuj) can also contain local (momentum 
independent) contribution which obviously vanishes in the limit of infinite dimensionality 
d — > oo and is not taken into account within the standard DMFT. Due to this fact there 
is no double counting of diagrams problem within our approach for the Hubbard model. 
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This question does not arise at all if we consider Ek(i&>) appearing due to some "additional" 
interaction. More important is that the assumed additive form of self-energy E(iu/) + Ek(icu) 
implicitly corresponds to neglect of possible interference of these local (DMFT) and non- 
local contributions. 

The self-consistency equations of our generalized DMFT+Ek approach are formulated 
as follows^: 

1. Start with some initial guess of local self-energy e.g. H(iu) = 0. 

2. Construct Ek(iu;) within some (approximate) scheme, taking into account interactions 
with collective modes or order parameter fluctuations which in general can depend on 
E(iu;) and \i. 

3. Calculate the local Green function 

GuM = ^11^ + ^- £ ( k ) _ E {iu) - E k (iw) ' (2) 

4. Define the "Weiss field" 

^'H^h + g^m. (3) 

5. Using some "impurity solver" to calculate the single-particle Green function for the 
effective Anderson impurity problem, defined by Grassmanian integral 

G d (r -t>) = 1-\ DctDc ia cUr)ci(r') exp(-S cff ) (4) 

^eS J 

with effective action for a fixed site ("impurity") % 

,/3 r /3 



rP rP rP 

Ses=- / dr\ j dT 2 Cia{ri)So l {Ti - r 2 )c+(r 2 ) + / drUn^(T)iii 
Jo Jo Jo 



(5) 



Z c s = j Dcf a Dci a exp(—S e ff), and (3 = T 1 . This step produces a new set of values 



G- d \iu). 



6. Define a new local self-energy 

V{iu) = Q Q - 1 {iu)-G- d \iu). (6) 
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7. Using this self-energy as "initial" one in step 1, continue the procedure until (and if) 
convergency is reached to obtain 



G u (iuj) = G d {iu). (7) 

Eventually, we get the desired Green function in the form of (Q), where Y,(iuj) and ^(iu) 
are those appearing at the end of our iteration procedure. 

III. CONSTRUCTION OF K DEPENDENT SELF ENERGY 

For the momentum dependent part of the single-particle self-energy we concentrate on 
the effects of scattering of electrons from collective short-range SDW-like antiferromagnetic 
spin (or CDW-like charge) fluctuations. To calculate Y>k(iuj) for an electron moving in 
the quenched random field of (static) Gaussian spin (or charge) fluctuations with domi- 
nant scattering momentum transfers from the vicinity of some characteristic vector Q ( "hot 
spots" model 2 ), we use a slightly generalized version of the recursion procedure proposed 
in Refs . 4 i 5 i 12 which takes into account all Feynman diagrams describing the scattering of 
electrons by this random field. This becomes possible due to a remarkable property of our 
simplified version of "hot spots" model that the contribution of an arbitrary diagram with 
intersecting interaction lines is actually equal to the contribution of some diagram of the 
same order without intersections of these line^^. Thus, in fact we can limit ourselves to 
consideration of only diagrams without intersecting interaction lines, taking the contribu- 
tion of diagrams with intersections into account with the help of additional combinatorial 
factors, which are attributed to "initial" vertices or just interaction lines^. As a result we 
obtain the following recursion relation (continuous fraction representation^ 2 -) for the desired 
self-energy: 

£ k (iu;) = £ n= i(iu;k) (8) 

with 

E ^ k ) = A2 —: — — r^r- v i v\ ■ (9) 

iuj + n — L{tuj) — e n {k) + inv n K — L n+ i(2Co>k) 
The quantity A characterizes the energy scale and n = C," 1 is the inverse correlation length 
of short range SDW (CDW) fluctuations, e n (k) = e(k + Q) and v n = \v£. + q\ + l^k+ol f° r 
odd n while e n (k) = e(k) and v n = \v£\ + for even n. The velocity projections v£ and 
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are determined by usual momentum derivatives of the "bare" electronic energy dispersion 
e(k). Finally, s(n) represents a combinatorial factor with 



s{n)=n (10) 

for the case of commensurate charge (CDW type) fluctuations with Q = (n/a, ir/a)^. For 
incommensurate CDW fluctuations^ one finds 

, s f 2±i for odd n 
s{n)=l 2 (11) 

| for even n. 

If we take into account the (Heisenberg) spin structure of interaction with spin fluctuations 
in "nearly antiferromagnetic Fermi-liquid" (spin-fermion (SF) model Ref.— ), the combina- 
torics of diagrams becomes more complicated. Spin-conserving scattering processes obeys 
commensurate combinatorics, while spin-flip scattering is described by diagrams of incom- 
mensurate type ("charged" random field in terms of Ref. 4 ). In this model the recursion 
relation for the single-particle Green function is again given by (JHJ), but the combinatorial 
factor s(n) now acquires the following form 4 -: 

for odd n 

(12) 

| for even n. 

Obviously, with this procedure we introduce an important length scale £ not present in 
standard DMFT. Physically this scale mimics the effect of short-range (SDW or CDW) 
correlations within fermionic "bath" surrounding the effective Anderson impurity. We expect 
that such a length-scale dependence will lead to a competition between local and non-local 
physics. 

An important aspect of the theory is that both parameters A and £ can in principle 
be calculated from the microscopic model at hand. For example, using the two-particle 
selfconsistent approach of Ref.-^ with the approximations introduced in Refs4»^, one can 
derive 1 ^ within the standard Hubbard model the following microscopic expression for A: 

A^ = -U 2 < n iT > + < nn > -2 < n iT n u > = 

U 2 < (rii T - nn) >= 



rr 



u2 <^>i < §1 >t (13) 



where we consider only scattering from antiferromagnetic spin fluctuations. Different local 
quantities here - spin fluctuation < Sf >, density n and double occupancy < n^n^ > - 
can easily be calculated within the standard DMFT— . We performed such calculations 11 for 
wide range of U and filling factors n using quantum Monte-Carlo (QMC) 15 . From these 
calculations we can see that the values of A lie in the interval of A = (0.5-i-2.0)£ and change 
rather smoothly with n and U. 

Microscopic expressions for the correlation length £ = /t _1 can also be derived within the 
two-particle self-consistent approach^. However, we expect those results for £ to be less 
reliable, because this approach is valid only for relatively small (or medium) values of U/t 
and for purely two-dimensional case (while real systems are quasi-two-dimensional). 

Thus, in the following we will consider both A and especially £ as some phenomenological 
parameters to be determined from experiments. This makes our approach somehow similar 
in the spirit to Landau approach to Fermi-liquids. 

Our construction can be further generalized to include other types of interactions. Thus 
scattering by random impurities with point - like potential V is easily taken into account in 
self - consistent Born approximation^. Then, in comparison with impurity free case, just 
we have a substitution (renormalization): 

e n -> e n - pV 2 ^ ImG(e n p) = e n r) e (14) 
p 

i, e =l-^^/mG( £nP ) (15) 
e n p 

If we do not perform fully self - consistent calculations of impurity self - energy, in the 
simplest approximation we just have: 

£n -> SnVe = £n + ^sigUEn (16) 

Ve=l + y L i (17) 
| 

where 7 = irpV 2 N (0) is the standard Born impurity scattering rate (AT (0) is the density 
of states of "free" electrons at the Fermi level). 
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IV. RESULTS AND DISCUSSION 



A. Computation details 

In the following, we discuss results for a standard one-band Hubbard model on a square 
lattice. With nearest (t) and next nearest (t') neighbour hopping integrals the dispersion 
reads 

£(k) = — 2t(cos k x a + cos k y a) — At' cos k x a cos k y a , (18) 

where a is the lattice constant. The correlations are introduced by a repulsive local two- 
particle interaction U. We choose as energy scale the nearest neighbour hopping integral t 
and as length scale the lattice constant a. 

For a square lattice the "bare" bandwidth is W = 8t. To study a strongly correlated 
metallic state obtained as doped Mott insulator we use U = AOt as value for the Coulomb 
interaction and a filling n = 0.8 (hole doping). The correlated metal in the case of W > U 
is considered for the case of U = At and filling factor n = 0.8 (hole doping). For A we 
have choosen rather typical values between A = O.lt and A = 2t (actually as approximate 
limiting values obtained from (fTT?J) via QMC calculations in Ref-^) and for the correlation 
length we considered mainly £ = 2a and £ = 10a (being motivated mainly by experimental 
data for cuprates^). 

The DMFT maps the lattice problem onto an effective, self-consistent impurity defined 
by Eqs. (0}-©- In our work we employed as "impurity solvers" two reliable numerically 
exact methods — quantum Monte-Carlo (QMC) 15 and numerical renormalization group 
(NRG) 1 ? 1 ^ . Calculations were done both for t' = and t'/t=-0.A (more or less typical for 
cuprates) at two different temperatures T = 0.088t and T = 0.356t (for NRG computations). 
QMC computations of double occupancies as functions of filling were done at temperatures 
T = O.lt and T = O.At. 

Below we present results only for most typical dependences and parameters, more details 
can be found in Ref.— . 
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B. Generalized DMFT+Ek approach: densities of states 

Let us start the discussion of results obtained within our generalized DMFT+Sk ap- 
proach with the densities of states (DOSs) for the case of small (relative to bandwidth) 
Coulomb interaction U = At with and without pseudogap fluctuations. As already discussed 
in the Introduction, the characteristic feature of the strongly correlated metallic state is the 
coexistence of lower and upper Hubbard bands split by the value of U with a quasiparti- 
cle peak at the Fermi level. Since at half-filling the bare DOS of the square lattice has a 
Van-Hove singularity at the Fermi level (t! = 0) or close to it (in case of t'/t = —0.4) one 
cannot treat a peak on the Fermi level simply as a quasiparticle peak. In fact, there are 
two contributions to this peak from (i) the quasiparticle peak appearing in strongly cor- 
related metals due to many-body effects and (ii) the smoothed Van-Hove singularity from 
the bare DOS. In Fig. ^ we show the corresponding DMFT(NRG) DOSs without pseudogap 
fluctuations as black lines for n = 0.8 for both bare dispersions t'/t = —0.4 (left panels) 
and for t' — (right panels) for two different temperatures T = 0.356t (middle panels) and 
T = 0.088t (upper and lower panels). The remaining curves in Fig. ^ represent results for 
the DOSs with non-local fluctuations switched on. For all sets of parameters one can see 
that the introduction of non-local fluctuations into the calculation leads to the formation of 
pseudogap within the quasiparticle peak. 

For n = 0.8 (Fig. ^) the picture of DOS is slightly asymmetric. The width of the pseudo- 
gap (the distance between peaks closest to Fermi level) appears to be of the order of ~ 2 A. 
We have checked that decreasing the value of A from 2t to t leads to a pseugogap that is cor- 
respondingly twice smaller and in addition more shallow. When one uses the combinatorial 
factors corresponding to the spin-fermion model fEq. (jl2|0 . the pseudogap becomes more 
pronounced than in the case of commensurate charge fluctuations (combinatorial factors of 
Eq. (II The influence of the correlation length £ is also as expected. Changing 
from = 0.1 to = 0.5, i.e. decreasing the range of the non-local fluctuations, slightly 
washes out the pseudogap. Also, increasing the temperature from T = 0.088t to T = 0.356t 
leads to a general broadening of the structures in the DOSs. Noteworthy is the fact that 
for t'/t = —0.4 and = 0.5 the pseudogap has almost disappeared for the temperatures 
studied here. Also very remarkable point is the similarity of the results obtained with the 
generalized DMFT+Ek approach with U — At (smaller than the bandwidth W) to those 
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obtained earlier without Hubbard-like Coulomb interactions^. 

Let us now consider the case of a doped Mott insulator. The model parameters are again 
taken as t'/t = —0.4 with filling factor of n = 0.8, but the Coulomb interaction strength is 
set to U — 401 Characteristic features of the DOS for such a strongly correlated metal are a 
strong separation of lower and upper Hubbard bands and a Fermi level crossing by the lower 
Hubbard band (for non-half-filled case). Without non-local fluctuations the quasi-particle 
peak is again formed at the Fermi level, but now the upper Hubbard band is far to the 
right and does not touch the quasiparticle peak (as it was for the case of small Coulomb 
interactions). 

Pseudogap appears close to the middle of quasiparticle peak. In addition we observe that 
the lower Hubbard band is slightly broadened by fluctuation effects. Qualitative behaviour 
of the pseudogap anomalies is similar to those described above for the case of U = 4t, e.g. 
a decrease of £ makes the pseudogap less pronounced, while reducing A from A = 2t to 
A = t narrows of the pseudogap and also makes it more shallow. Note that for the doped 
Mott-insulator the pseudogap is remarkably more pronounced for the SDW-like fluctuations 
than for CDW-like fluctuations. 

There are, however, quite clear differences to the case of U — At. For example, the width 
of the pseudogap appears to be much smaller than 2A, which we attribute to the fact that 
the quasiparticle peak itself is actually rather narrow in the case of doped Mott insulator. 

Random impurity scattering, in general case leads to the filling of the pseudogap with 
the growth of impurity scattering rate both for correlated metal and doped Mott insulator. 
As a typical example, in Fig. 01 we show results of our calculations for the case of doped 
Mott insulator. These were obtained via non self-consistent procedure (using (fTfij). (JHJ), as 
full self-consistent procedure leads only to rather insignificant quantitative changes. 



C. Generalized DMFT+Sk approach: spectral functions A(uj,\s) 

In the previous subsection we discussed the densities of states obtained self-consistently 
by the DMFT+Sk approach. Once we get a self-consistent solution of the DMFT+Ek 
equations with non-local fluctuations we can, of course, also compute the spectral functions 
A(u,k) 

A(u,k) = —lm — — ^ — — — , (19) 
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where self-energy £(o>) and chemical potential n are calculated self-consistently as described 
in Sec. |H] To plot A(u>,\t) we choose k-points along the "bare" Fermi surfaces for different 
types of lattice spectra and fillings. In Fig. |U one can see corresponding shapes of these 
"bare" Fermi surfaces (presented are only 1 / 8-th parts of the Fermi surfaces within the first 
quadrant of the Brillouin zone). 

In the following we concentrate mainly on the case U = 4t and filling n = 0.8 (Fermi 
surface of Fig. Ufa)). The corresponding spectral functions A(u,k) are depicted in Fig. |SJ 
When t'/t = —0.4 (upper row), the spectral function close to the Brillouin zone diagonal 
(point B) has the typical Fermi-liquid behaviour, consisting of a rather sharp peak close to 
the Fermi level. In the case of SDW-like fluctuations this peak is shifted down in energy 
by about — 0.5t (left upper corner). In the vicinity of the "hot-spot" the shape of A(u;,k) 
is completely modified. Now A(u, k) becomes double-peaked and non-Fermi-liquid-like. 
Directly at the "hot-spot", A(u,\i) for SDW-like fluctuations has two equally intensive 
peaks situated symmetrically around the Fermi level and split from each other by ~ 1.5 A 
Refs.-^. For commensurate CDW-like fluctuations the spectral function in the "hot-spot" 
region has one broad peak centered at the Fermi level with width ~ A. Such a merging 
of the two peaks at the "hot-spot" for commensurate fluctuations was previously observed 
in Ref.~. However close to point A this type of fluctuations also produces a double-peak 
structure in the spectral function. 

In the lower panel of Fig. El we show spectral functions hole doping (n = 0.8) and the 
case of t' = (Fermi surface from Fig. Efb)). Since the Fermi surface now is everywhere 
close to the antiferromagnetic zone boundary, the pseudogap anomalies are rather strong 
and almost non-dispersive along the Fermi surface. 

For the case of a doped Mott insulator (U = 40i, n = 0.8), the spectral functions obtained 
by the DMFT+E k approach are presented in Fig. |BJ Qualitatively, the shapes of these 
spectral functions are similar to those shown on Fig. 03 As was pointed out above, the strong 
Coulomb correlations lead to a narrowing of the quasiparticle peak and a corresponding 
decrease of the pseudogap width. One should also note that in contrast to U = At the 
spectral functions are now less intensive, because part of the spectral weight is transferred 
to the upper Hubbard band. 

Using another quite common choice of k-points we can compute A(u, k) along high- 
symmetry directions in the first Brillouin zone: T(0, 0)— X(7r, 0)— M(n, 7r)— T(0, 0). The spectral 
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functions for these k-points are shown in Fig. [7|for the case of SDW-like fluctuations and 
U = At. For all sets of parameters one can see a characteristic double - peak pseudogap 
structure close to the X point. In the middle of M — T direction (so called "nodal" point) 
one can see the reminiscence of AFM gap which has its biggest value here in case of perfect 
antiferromagnetic ordering. Also in the nodal point "kink" -like behavior is observed caused 
by interactions between correlated electrons with pseudo-gap fluctuations. A change of the 
filling leads mainly to a rigid shift of spectral functions with respect to the Fermi level. For 
the case of U — 40t spectral densities demonstrate rather similar behavior—. 

With the spectral functions we are now in a position to calculate angle resolved photoe- 
mission spectra (ARPES), which is the most direct experimental way to observe pseudogap 
in real compounds. For that purpose, we only need to multiply our results for the spec- 
tral functions with the Fermi function at appropriate temperature. Typical example of the 
resulting DMFT+Ek ARPES spectra are presented in Fig. |HJ One should note that for 
t'/t = —0.4 (upper panel of Fig. (HJ) as k goes from point "A" to point "B" the peak situated 
slightly below the Fermi level changes its position and moves down in energy. Simultane- 
ously it becomes more broad and less intensive. The dotted line guides the motion of the 
peak maximum. Such behaviour of the peak in the ARPES is rather reminiscent of those 
observed experimentally in underdoped cuprates^ 4 * 1 *- 

D. "Destruction" of the Fermi surface 

Within the standard DMFT approach Fermi surface is not renormalized by interactions 
and just coincides with that of the "bare" quasiparticles^. However, in the case of nontrivial 
momentum dependence of electron self - energy, important renormalization of the Fermi 
surface appears due to pseudogap formation 4 . There are a number of ways to define Fermi 
surface in strongly correlated system with pseudogap fluctuations. In the following we are 
using intensity plots (within the Brillouin zone) of the spectral density (JT§j) taken at u = 0. 
These are readily measured by ARPES and appropriate peak positions define the Fermi 
surface in the usual Fermi liquid case. 

Our results 20 are shown in Fig. H3 for the case of correlated metal with U — At and in 
Fig. El for the doped Mott insulator (U = 40i) (in both cases we assume spin - fermion 
combinatorics). The qualitative behavior observed in Fig. |H1 clearly demonstrates the "de- 
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struction" of the well defined Fermi surface in the strongly correlated metal with the growth 
of the pseudogap amplitude A. Quite similar behavior was first observed in pioneering paper 
by Norman et al.— and in numerous later ARPES experiments. It is seen, that "destruc- 
tion" of the Fermi surface starts in the vicinity of "hot spots" for small values of A, but 
almost immediately it disappears in the whole antinodal region of the Brillouin zone, while 
only "Fermi arcs" remain in the nodal region very close to the "bare" Fermi surface. These 
results give a natural explanation of the observed behavior and also of the fact that the 
existence of "hot spots" regions was observed only in some rare cases^. 

For the case of doped Mott insulator shown in Fig. El we see that the Fermi surface 
is rather poorly defined for all values of A, as the spectral density profiles are much more 
"blurred" than in the case of smaller values of U, reflecting important role of correlations. 

It is interesting to note that from Figs.OHEIlit is clearly seen that rather natural definition 
of the Fermi surface as defined by the solution of the equation 

w-e(k) + fx-ReE(u) - ReE k (u) = (20) 

for u = 0, used e.g. in Ref.» 4 is inadequate for strongly correlated systems with finite U and 
nonlocal interactions (pseudogap fluctuations). 

V. CONCLUSION 

To summarize, we propose a generalized DMFT+Ek approach, which is meant to take 
into account the important effects due to non-local correlations in a systematic, but to 
some extent phenomenological fashion. The main idea of this extension is to stay within 
a usual effective Anderson impurity analogy, and introduce length scale dependence due to 
non-local correlation via the effective medium ("bath") appearing in the standard DMFT. 
This becomes possible by incorporating scattering processes of fermions in the "bath" from 
non-local collective SDW-like antiferromagnetic spin (or CDW-like charge) fluctuations. 
Such a generalization of the DMFT allows one to overcome the well-known shortcoming of 
k-independence of self-energy of the standard DMFT. It in turn opens the possibility to 
access the physics of low-dimensional strongly correlated systems, where different types of 
spatial fluctuations (e.g. of some order parameter), become important. However, we must 
stress that our procedure in no way introduces any kind of systematic 1 / rf-expansion, being 
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only a qualitative method to include length scale into DMFT. 

In our present study we addressed the problem of pseudogap formation in the strongly 
correlated metallic state. We showed evidence that the pseudogap appears at the Fermi level 
within the quasiparticle peak, introducing a new small energy scale of the order of A in the 
DOSs and spectral functions A(u},k) and significant renormalization of the Fermi surface. 

Let us stress, that our generalization of DMFT leads to non-trivial and in our opinion 
physically sensible k-dependence of spectral functions. Similar results were obtained in 
recent years using the cluster mean-field theories^. The major advantage of our approach 
over these theories is, that we stay in an effective single-impurity picture. This means that 
our approach is computationally much less expensive and therefore also easily generalizable 
for the account of additional interactions . 
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FIG. 1: Comparison of DOSs obtained from DMFT(NRG)+£k calculations for different combi- 
natorical factors (SF — spin-fermion model, commensurate), inverse correlation lengths in 
units of the lattice constant, temperatures (T) and pseudogap potential A = 2t. Left column 
corresponds to t'/t = —0.4, right column to t' = 0. In all graphs the Coulomb interaction is U = 4i 
and n = 0.8. The Fermi level corresponds to zero. 
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FIG. 2: Comparison of DOSs obtained from DMFT(NRG)+S k calculations for t'/t = -0.4, T = 
0.088t, U = AOt and filling n = 0.8. 
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FIG. 3: Comparison of DOSs obtained from DMFT (NRG) + S k calculations for t'/t = -0.4, T = 
0.088i, U = iOt, n = 0.8 and for impurity scattering rate 7 = 1.2t. 




FIG. 4: 1/8-th of the bare Fermi surfaces for the different occupancies n and combinations (t,t') 
used for the calculation of spectral functions A(\t, uj). The diagonal line corresponds to the anti- 
ferromagnetic Brillouin zone boundary at half-filling for a square lattice with nearest-neighbours 
hopping only. The full circle marks the so-called "hot-spot" . 
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FIG. 5: Spectral functions A(k,uj) obtained from the DMFT(NRG)+$]i c calculations along the 
directions shown in Fig. |1J Model parameters were chosen as U = At, n = 0.8, A = 2t, £ -1 = 0.1 
and temperature T = 0.088i. The "hot-spot" k-point is marked as fat dashed line. The Fermi 
level corresponds to zero. 
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FIG. 6: Spectral functions A(k,u) obtained from the DMFT(NRG)+S k calculations for U = AOt, 
other parameters as in Fig. |3 



For all graphs SF , 
M, U=4t, T=0.088t, n=0.8, 5~'=0.1 
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FIG. 7: Spectral functions A(k,u)) obtained from the DMFT(NRG)+$]k calculations along high- 
symmetry directions of first Brilloin zone r(0, 0)— X(7r, 0)— M(ir, tt)— T(0, 0), SF combinatorics (left 
row) and commensurate combinatorics (right column). Other parameters are U = At, n = 0.8, 
A = 2t, = 0.1 and temperature T = 0.088t. The Fermi level corresponds to zero. 




FIG. 8: ARPES spectra obtained from the DMFT(NRG)+S k calculations for U = At and n = 0.8 
along the lines in the first BZ as depicted by Fig. 0J all other parameters as in Fig. [SJ 




FIG. 9: Fermi surfaces obtained from the DMFT(NRG)+S k calculations for U = At and n = 0.8. 
Shown are intensity plots of spectral density (|19|) for u = and ^ = 10a. (a) - A = 0.2t; (b) - 
A = OAt; (c) - A = t; (d) - A = 2t. Dashed line denotes "bare" Fermi surface. Black lines show 
the solution of Eq. ((20)) . 
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FIG. 10: "Destruction" of the Fermi surface obtained from the DMFT(NRG)+£k calculations for 
U = 40t and n = 0.8. Other parameters and notations are the same as in Fig. [§J 



